Kink nucleation in two-dimensional Frenkel-Kontorova model 
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A computer simulation of thermofluctuation nucleation of kinks on dislocations and their dynamics 
is carried out in the framework of two-dimensional (2D) Frenkel-Kontorova model. It is shown that 
at relatively low temperatures and applied stresses the kinks can appear as a result of developing 
instability of phonon modes localized in the vicinity of the dislocation. The transition from this 
mechanism to the ordinary thermofluctuation kink nucleation with temperature increase can reveal 
itself in the peculiarities of yield stress temperature dependence. 
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I. INTRODUCTION 

Due to the translational symmetry of crystals the en- 
ergy of a dislocation depends periodically on its posi- 
tion in the lattice (Peierls relief )Ero. This concept having 
been formulated more than 50 years agctl proved to be 
very fruitful for the dislocation theory, in particular, for 
elucidations of factors determining the dislocation mobil- 
ity (see, e.g., the reviewscttl). At present it is commonly 
accepted that at finite temperatures the dislocations in 
crystals glide in the Peierls relief due to nucleation of 
kinks and their propagation along the dislocation lines. 
This mechanism determines .-the dislocation mobility in 
covalent crystals (like silicon)a but can also be essential in 
the bec metals and some intermetallic compounds (TiAl, 
NiAl)Q. 

According to the conventional point of view the veloc- 
ity of dislocations in the Peierls relief is determined by 
the balance of the processes of thermofluctuation nucle- 
ation of kinks and annihilation due to mutual recombi- 
nation of kiftks and antikinks or their disappearance at 
the defectsO'Q. HowevgiL recently the experimental data 
have become availableOu'El which cannot be described ad- 
equately by this simplified scheme. The resultsa of nu- 
merical simulation of dislocation migration in the 2D 
Frenkel-Kontorova model also have demonstrated that 
the traditional views should be revised. As it was shown 
in RefB the kinks on the dislocations can behave like soli- 
tons and, when collising, pass through each other without 
annihilation in spite of thermal fluctuations and damp- 
ing. As a result a much higher kink density than that 
predicted in terms of the traditional approach should be 
expected. It seems interesting to consider the question on 
the contribution of dynamical effects as well to the dou- 
ble kink nucleation since it is the process of kink-antikink 
pair nucleation that determines the temperature depen- 
dence of the dislocation mobilityu. 

In the present work the process of kink nucleation is 
thoroughly studied by means of modeling the dislocation 
movement in the 2D Frenkel-Kontorova model at the fi- 
nite temperatures. We demonstrate that, in addition to 



the well known mechanism of thcrmoflucutation nucle- 
ation of kinks, they may also appear as a result of devel- 
opment of instability of lattice vibrations localized near 
the dislocation axis. 



II. 



FORMULATION OF THE MODEL 



Similar to Ref.EI, we study the dislocation dynamics 
in terms of the 2D Frenkel-Kontorova (FK) model rfik-fi- 
nite temperatures. In the generalized 2D FK modelaO a 
layer of atoms, set into the periodic potential relief and 
interacting with each other by means of elastic forces, is 
considered. The potential energy of this system has the 
following form 
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where K is the stiffness of interatomic bonds, (n, m) des- 
ignates the sum over all the pairs of the nearest neighbors, 
gi = JL(1; 0), g2i3 = ^(|; ±% are the reciprocal lat- 
tice vectors of minimum length. 

In order to consider the damping and thermal fluctu- 
ations we introduce the contact with a thermostat using 
the method of Langevin equations of motionO 
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where u„ is the vector of displacement of the n-th atom 
from its equilibrium position, 7 is the friction coeffi- 
cient, i n is the external force assumed to be equal for 
all the atoms, (, n i(t) is the random Gaussian variable 
having the properties (£„»(*)) = 0, (£ni(*)C'i' (*')) = 
^T^nn' — * )' (* = x tV)- The brackets indicate av- 
eraging over the realization of random process £, n {t)- A 
method of solving the set ofrStochastic differential equa- 
tions (||) is described in Ref.B. 

For the numerical experiment the 2D hexagonal lattice 
with 40x40 atoms with periodic boundary conditions was 
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choosen. At the initial moment the displacements in the 
crystallite were specified in accordance with the known 
solution for dislocations in the ID continual FK modelH: 

1 , 35, , 1 [K 

u x = 0, u y — 1 H — arctanexp(— — J, A = —\ — 

7T A 7T V P 

(3) 

Then the set of equations was integrated numerically for 
a rather long time (t ~ 10 2 ) at f„ = and the choosen 
temperature T to find an equilibrium configuration of 
the dislocation. After that an external force f = (0, f y ), 
amounting to 0.75-^0.85 of limiting value fp correspond- 
ing to the Peierls stress, was applied. 

The most important parameter in the model ([!]), (|^) is 
the ratio of the height of potential relief P to the stiffness 
K . We put K = 1 (which determines the energy units) 
and choose P « 0.1 when, for small displacement u, the 
forces of atomic interaction with the neighbors and with 
the substrate layer are comparable. This ratio of param- 
eters leads to the dislocation with A ~ 1, which is typical 
for most of metals (seeE3ll3). The temperature T varied 
over a wide range from 10 -3 to 10~ 2 i.e. it was of or- 
der of (0.1 -r 0.01)P, and the friction coefficient 7 varied 
over the range (2 -=-5) • 10~ 2 , which corresponds to typical 
values for metals at room temperatures (see, e.g., Ref.El). 

III. COMPUTATIONAL RESULTS 

Fig.|l| shows a sequence of dislocation states in the 
2D lattice, demonstrating thermofluctuation nucleation 
of kink-antikink pairs. It is seen that their appearance 
is preceded by the stage of small amplitude oscillations 
of the dislocation segment. At a certain time one of 
these oscillation modes begins to increase, its amplitude 
reaches the maximum value equal to the distance be- 
tween the atom rows, and a kink-antikink pair is formed. 
At subsequent moments of time the kinks move along 
the dislocation line, reach the crystallite boundaries and, 
because of periodic boundary conditions, appear on its 
opposite edges. As it was demonstrated in Ref.cl, under 
the conditions when the termofluctuation kink nucleation 
is possible, the kinks demonstrate a soliton-like character 
of motion and, when meeting, pass through each other 
without annihilation. 

In order to understand the mechanism of the pair kink 
nucleation, let us consider in more detail the picture of 
atom displacements. An important feature, accompany- 
ing the kink nucleation at low values of external force f 
(or of the temperature), is the localization of long- wave 
lattice oscillations in the vicinity of the dislocation which 
manifests itself as a correlated motion of atoms in the 
rows nearest to the dislocation axis. This feature is illus- 
trated by the dependence of the position of center of mass 
of atomic rows on time (Fig.|]a). Curves 1,2 correspond 
to the y-component (parallel to the Burgers vector) of the 
average displacement < u y > for the atomic rows with 
n = 20 and 19 (see Fig.|l|) nearest to the dislocation axis, 



curves 3,4 for the atomic rows with n = 21 and 18 follow- 
ing the nearest ones. One can see from FigJ|a that the 
atom row n — 20 demonstrates rather regular displace- 
ment oscillations < u y > around the average value with a 
period of tens of characteristic phonon times. The ampli- 
tude of these oscillations decreases rapidly with the dis- 
tance from the dislocation axis, and for n = 18 it is prac- 
tically comparable with the thermal noise level (curve 4, 
Fig. ||a). Note that the amplitude of the ^-component 
of average displacements < u y > (perpendicularly to the 
Burgers vector) for rows n = 20 and 19 is more than by 
an order of magnitude smaller than < u y > . This reduces 
effectively the dimensions of the problem. The possibil- 
ity of localization of lattice oscillations on the disloca- 
tion was predicted earliercl in terms of a simple continual 
model. Here this phenomenon was found from the nu- 
merical experiment on the discrete lattice with allowance 
for the temperature and damping. 

The oscillation stage ends with the nucleation of kink 
pair whose propagation along the dislocation leads to 
monotonous growth in < u y > dependence on t (curve 1, 
Fig||) . Thus the analysis of atom displacements reveals 
that the kink pair nucleation is due to the loss of stabil- 
ity of atom row motion in the effective periodic potential 
and the onset of an inhomogeneous state similar to the 
crowdion-anticrowdion pair (see Fig.^). A similar._Dhc- 
nomenon for the ID FK model was discussed in Refo. It 
is this configuration in the row of atoms that corresponds 
to appearance of the kink-antikink pair on the disloca- 
tion line (Fig.|l|). Thus, the 2D calculations demonstrate 
clearly that there is a region of parameters where the sce- 
nario of kink pair nucleation on a dislocation is due to 
the oscillatory instability near the dislocation axis rather 
than to the classical pattern of thermoflucutation forma- 
tion and growth of the nucleus. In addition, the kinks 
on the dislocation in the 2D FK model arc found to be 
similar to crowdions in the ID FK model. Note that 
the results obtained not only support this commonly ac- 
cepted analogy but clarify and extend it. 

The observed picture of the kink pair nucleation has 
an uncommon non-monotonous dependence in parame- 
ters T, f (Fig.^J): with the rise on temperature and/or 
the magnitude of applied force the probability of kink 
nucleation first drastically reduces (in region II the kinks 
do not practically appear during the simulation) and with 
further increase in the parameters are easily formed again 
(region III). The diagram shown in Fig. [| is constructed 
for fixed parameters P and 7; the increase in the damp- 
ing results in the reduction of the region I. The curves in 
Fig.^b illustrate the system behavior in three different re- 
gions. Curve 1 corresponds to the kink nucleation at low 
temperature and force applied, described above in detail 
(region I). Curve 2 corresponds to the intermediate region 
II where the amplitude of oscillations < u y > decreases 
rapidly with time to a value comparable with the level of 
thermal noise. Curve 3 is typical of high temperatures or 
external forces (region III) when the kink nucleation is 
again observed. However the "high-temperature" kinks 
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differ essentialiy from the "low-temperature" ones; they 
are much wider so that displacements change slowly at 
a distance to some lattice parameters (Fig |3j3) . Their 
nucleation is more similar to the classical pattern of fluc- 
tuation nucleation and growth of the nucleate, without 
the stage of preliminary oscillations of the dislocation line 
(Fig.|l|) which is typical of region I. 

IV. DISCUSSION 

The computer simulations in the 2D FK model de- 
scribed above demonstrate a new mechanism of kink nu- 
cleation on a dislocation as a result of development of 
the instability in the phonon subsystem. Along with 
the phonon modes spreading within the whole crystallite 
(bulk phonons) there are also modes localized near the 
dislocation. In terms of continuum model the possibility 
of existence of gpch localized oscillations was discussed 
earlier (see Reftl). These oscillations are clearly seen in 
Fig.|. 

At not too high temperatures and forces (region I, 
Fig.^) the motion of atoms in the rows nearest to the 
dislocation line appears to be correlated, which manifests 
itself in regular oscillations < u y > (Fig|2|a) . The change 
in the behavior of atom motion in the transition to the re- 
gion II of "intermediate" temperatures and forces (curve 
2, FigJ|) is accounted for by the increased damping of 
localized phonons and loss of coherency, first in the dis- 
placements of neighbor atom rows and then inside them. 
As a result the kink nucleation turns out to be suppressed 
(note that the growth of damping leads to the same con- 
sequences). With further increase in the temperature 
(curve 3, Fig.||) the kinks again begin to nucleate but in 
this case it is a result of development of the instability of 
"common" bulk phonons, which corresponds to the tradi- 
tional mechanism of thermoflucutation nucleation of the 
double kinktl. The transition from region I to region III 



is accompanied by the change in the kink structure, from 
narrow kinks (Fig |3]a) to the wide ones (Fig.|]b). This in- 
dicates the decrease in the secondary Peierls relieia. The 
latter is that shows the loss of coherence in the atom 
motion in transition to "high-temperature" region III. 

Let us discuss the significance of the results obtained 
for understanding the regularities of dislocation motion 
in the Peierls relief. This motion, generally speaking, 
occurs due to a complicated collective process with the 
participation of many kinks. We confined ourselves here 
to the consideration of the initial stage: nucleation of 
single kink-antikink pairs. A new phenomenon of fast 
kink nucleation (for the characteristics times of order of 
10-100 inverse phonon frequencies) has been found. In 
reality, apart from this process a slower thermoactivated 
kink nucleation process takes place. In the present work 
we did not study the temperature dependence of the kink 
nucleation rate. One can assume, however, that for the 
fast processes considered here, this dependence is weaker 
than for slow termoactivated processes. Therefore one 
can expect that the transitions from the low temperature 
region to the intermediate temperature one (suppression 
of the process of kink nucleation) and from the intermedi- 
ate to high temperature region (restarting of this process) 
may correspond to inflection points of the temperature 
dependence of yield stress which have been actually ob- 
served in sorae_metals and alloys with a sufficiently high 
Peierls reliefoEj. However, this question needs further 
investigations. 
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FIG. 1: Configuration of the dislocation line (determined by the condition \u y \ = 0.5b) at different times for P = 0.1, T = 0.01, 
7 = 0.02, / = —0.061. The rows closest to the dislocation line are n = 19, 20. 
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FIG. 2: Time dependence of the center of mass dispacement of atomic chains for the rows nearest to the dislocation line with 
n = 20, 21, 19, 18 (curves 1-4) at P = 0.1, 7 = 0.02 T = 0.01, / = -0.061 (a) and for atom row with n = 20 P = 0.1, 7 = 0.02, 
T = 0.01, / = -0.061 (curve 1), T = 0.08, / = -0.07 (curve 2), T = 0.01, / = -0.07 (curve 3) (b). Solid line is the result of 
smearing of the computational data to eliminate the thermal noise. 
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FIG. 3: Crystallite fragments showing the structure of "low-temperature" (T = 0.003, f y = —0.07, a) and "high-temperature" 
(T = 0.01, f y — —0.068, b) kinks. The regions of the compression and expansion corresponds to crowdion and anticrowdion, 
correspondingly. 
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FIG. 4: Regions of the parameter values where the kink-antikink pair nucleation occurs from localized oscillations (I), kinks 
do not nucleate (II), and they nucleate by the standard thermofiuctation mechanism (III) at P — 0.1 7 = 0.02. The circles 
show the parameter values at which the change in the mechanism was observed. 



